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Motivation: Sign problem, Silver Blaze problem 
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Finite-density quantum chromodynamics (QCD) 


QCD 
Fundamental theory for quarks and gluons 


Neutron star 


e Cold and dense nuclear matter 


o 2Msun neutron star (2010) 


Neutron star merger 
o Gravitational-wave observations (2016) (image from NASA) 


Reliable theoretical approach to equation of state must be developed! 


Zaco(T,p)= [DA Det(P(A, pa) + m) exp C Svu(). 


quark gluon 


Sign problem: Det(D(A, uq) + m) Z 0 at uq 40. 
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Sign problem in finite-density QCD 


QCD & |QCD| 


ZQCD = [ra (det JD) e SYM | Ziacp| = In det P| eM, 


If these two were sufficiently similar, we have no practical problems. 
However, it was observed in lattice QCD simulation that at T = 0 


(e.g., Barbour et. al. (PRD 56 (1998) 7063)) 


nB QCD 


E > 
mM /2 -mn/3 CS 


Yuya Tanizaki (RIKEN, RBRC) Lefschetz-thimble path integral 


Jun 17, 2016 € Tokyo 5/38 


Baryon Silver Blaze problem 


The curious incident of the dog in the night-time (Holmes, Silver Blaze). 
Problem: It is of great importance for finite-density lattice QCD 
to understand ng = 0 for fg < My /3. (Cohen, PRL 91 (2003) 222001) 


Current situation: For u, < m„/2, the problem is almost solved: 
Quark det. is (Adams, PRD 70 (2004) 045002) 


Det(D(A, el +m) | 
Det(D(A,0) tm) — 


I sep äi, — Aa); 


Re(A4)<uq 


and ess-min4(Re(A4)) = Mx/2 (Nagata et. al., PTEP 2012 014103). 


For fg > Mz/2, no one knows how to understand this. 
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Method: Path integral on Lefschetz thimbles 
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Sign problem of path integrals 


Consider the path integral: 


ja EEGEN 


e S|x] is real = No sign problem. Monte Carlo works. 
e S|x] is complex = Sign problem appears! 


If S[x] € C, eom S'[x] = 0 may have no real solutions z(t) € R. 


Idea: Complexify x(t) € C! 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = | evi E + 2 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a — 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
d 3 
Ai(a) = > No [ Pi E taz) 


ng: intersection number of the steepest ascent contour Ke and R. 


"T n 


Figure : Lefschetz thimbles J and duals K 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent" cycles J,: (classical eom S'(2,) = 0) 


/ dx e- 560 = Y (Ko, y) d”z e80, 


o 


Js are called Lefschetz thimbles, and Im[S] is constant on it: 


(Ko, 


I= {0m 0-2), P 


IR): intersection numbers of duals K, and | 


(Ko = {2(0)|2(00) = zo 1). 
[Witten, arXiv:1001.2933, 1009.6032] 
[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
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Previous studies: MC with one-thimble ansatz 
In previous studies, one has been used one-thimble ansatz: 


i d"x eS) sf d^z e- 90. 
m Jo 


[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
Relativistic Bose gas: 


Thimble 2-(b) +8: " T 
3 © Langen" > 
25 
$ 2 
2 
3 
3 
5 15 
al 
D 
E 
v 
# E 
0.5 
z 
o x y , , f f 
0.9 1 44 1.2 1.3 1.4 1.5 1.6 
u (chemical potential in lattice unit) Fujii et al (JHEP 1310) 
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Previous studies: MC with one-thimble ansatz 
Chiral Random Matrix: 


d DUX dog Malen) 


N 
icosh(u) X + sinh(u)Y f 
x [det ee sinh(p1) Y t = m — ) 


20 


nao 


CRMT with 1-thimble ansatz with Nr = 2, u/VN = 2. (Di Renzo, Eruzzi, GE 015)) 
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Analysis: Semi-classical analysis of the one-site Hubbard model 
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One-site Fermi Hubbard model 
One-site Hubbard model: 


A = Uff, — u(fu + f). 


Fock state gives the number density immediately: 


Bai 


(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258.) 
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Analysis 


Path integral for one-site model 


The path-integral expression for the one-site Hubbard model: : 


B d A 2. UD 
a=a > | ay (14 ef Ge)? 9807/20, 
27U R 2 Ub cia addio D 


Fermion Det 


Integrand has complex phases causing the sign problem. 


( is an auxiliary field for the number density: 
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Sign problem and fermion determinant 


One-site Hubbard model: 
U gi -&(-U/2- u) oib)? 
Det |0, — jm enit = (l+te mete)”, 


Quark determinant in QCD: 


Det kaaft, +m) — u] = N (A) I] (14-678 6;7u-i6)) (]4e-Blertutien)), 


ej>0 
where the spectrum of ya (1D, + m) is 
Atıın) = Ej (A) = id; (A) + (2n + Du. 


Minimal value of e(4) = m; /2. 
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Silver Blaze problem for y < —U/2, u < m7/2 


One-site Hubbard model: As BU > 1 and —U/2 — u > 0, 


U ; 
Det 2, = (» + E + e) = (1 + EURER)" ~]. 


The sign problem almost disappears, so that J, ~ IR. 
Finite-density QCD: As 6 — oo and u < m,/2, 


Det ud, +m) — ul 
Det l4 D 4 F m)| 


The sign problem disappears by the reweighting method. 
=> Lefschetz thimbles ~ Original integration regions 


> 1 
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Analysis 


Flows at u/U < —0.5 (and u/U > 1/5) 


WH 
0 
Ge 
^ d 
| 
ies ll 
ep 
E 
u 
+ 
bd 
8 3 
2 
LL. 
© 
ui Pun 
> = 
.20 
EE, 


| dz e- 909, 
Se 


Number density: n, = 0 for u/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 


19 / 38 


Lefschetz-thimble path integral Jun 17, 2016 € Tokyo 


Yuya Tanizaki (RIKEN, RBRC) 


Silver Blaze problem for u > —U/2, u > m,/2 


One-site Hubbard model: At each real config., the magnitude is 
exponentially large: 


Det H — (» + > + e) =0 kk 


This large contributions must be canceled exactly in order for n = 0. 
Finite density QCD: The situation is almost the same, since 


Det(D(A, uj) tm) _ 
Det(D(A,0)+m) ` 


I ep gin — AA), 


Re(AA)<Hq 


but ng = 0 for u, S My /3. 
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Analysis 


Flows at —0.5 < u/U < 1.5 
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Flow at u/U — 0 


Figure : 


gU 1. 


This value is far away from n = Im (z)/U=0, 1, or 2. 


mn 


Complex saddle points lie on Im(z,,)/U 
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Curious incident of n in one-site Hubbard model 


We have a big difference bet. the exact result and naive expectation: 


05 0 0.5 1 15 ajo 


This is similar to what happens for QCD and |QCD|. 
p/U = -0.5 & u, = Mx /2. 
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Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


WU x 1/2 


U 
Ea cR By +2rmT. 


At these solutions, the classical actions become 


SE (B +3)’, 
Re (Sm — So) ~ 2m m?, 
Im S. c 2mm (5 4- 1). 
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Semiclassical partition function 


Using complex classical solutions Zm, let us calculate 


oo 
La := d ee, 


M=—-09 


This expression is valid for —1/2 S u/U $ 3/2. 
This is calculable using the elliptic theta function: 


= 1 
-S [149 9 Heo \ 20? pl 
e + d cos 27m (5 + 2 e 


Za = 
m=1 
” H 1 942 
= ep |r| ES |,6 2" IER 
U 2 
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Number density & Lee-Yang zeros 
Lee-Yang zeros of Za: 


e e 
e e 
a e Re 
e e 


Semiclassical study gives the correct transition! 


" 2 (1« u/U « 3/2), 
Nel = 39, Za + 1 (0 « p/U « 1), 
Pap 0 (-1/2<p/U <0). 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Important interference among multiple thimbles 


Let us consider a “phase-quenched” multi-thimble approximation: 


Za Aker Mriéale 0727). 


o Lee-Yang zeros cannot appear at u/U =0, 1. 
e One-thimble, or “phase-quenched”, result: n ~ u/U + 1/2. 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Particle Productions: Dumulu, Dunne, PRL 104 250402) 
(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schäfer, Ünsal, PRL 115 041601) 
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Numerical results 


Results for BU = 30: (1, 3, 5-thimble approx.: Jo, Jo U Ji, and Jo U J41 U Ja ) 


2.0 
— exact 
— 1-thimble 


15 --- 3-thimble 


i 
1 
" 
11 
[E 
14 
rn 
14 
UI 
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UI 
ry 
ry 
ry 
ry 
1 
1 


Necessary number of Lefschetz thimbles ~ OU (27). 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002, arXiv:1509.07146[hep-th]) 
Yuya Tanizaki (RIKEN, RBRC) 
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Bonus: Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz (0) _ 8S 
0: Stochastic time, n: Random force (n(0)n(0')),, = 20(0 — 0"). 
Ito calculus shows that 


x (Olen(0)) = A(O" (z,(8))), — (O'(Zn(8))S"(2n(9)))n- 


If the I.h.s becomes zero as 0 > oo, this is the Dyson-Schwinger eq. 
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Complex Langevin study of one-site Hubbard model 


BU =30 /U=0 


“AS  -10  -5 0 5 10 15 
BRe(z) 
Figure : CL distribution at u/U = 0 


It looks quite similar to Lefschetz thimbles, however the ensemble 
average cannot reproduce their complex phases. 
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Complex Langevin simulation 
One-site Fermi Hubbard model: 


2.5 
--:- exact 

2.0} HH w/o phase 
‚91 w/ phase 


-0.5 0.0 0.5 1.0 Lë 
p/U 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 


Consequence 


The knowledge of LT helps us to understand the longstanding failures 
of CL method. 
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Review 


Review: Practical algorithm for simulating multiple Lefschetz thimbles 
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Review 


Possible concerns for practical applications 


In this study, we emphasized the importance of interference among 
Lefschetz thimbles. 


Then, the following becomes important: 
e Find all contributing complex saddle points, 
e Construct Lefschetz thimbles around those saddle points, 
e Evaluate integration on each Lefschetz thimbles, and 
e Sum up those results. 
We need some machinery to do them automatically. 
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Review 


Neat idea for multiple thimble simulation 


d oS 
Deform the original cycle IR” by the gradient flow, E — (5) 


(Alexandru, et. al. JHEP (2016)) 


Im(z) 
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Review 


Neat idea for multiple thimble simulation 


d oS 
Deform the original cycle IR” by the gradient flow, E — (5) 


(Alexandru, et. al. JHEP (2016)) 


Im(z) 
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Review 


Formulation 
Let us fix a flow time T', and define 


Are far x) € C” ie = (5) z(0;::z) =a €] dë 


A 


By construction, z(T;-) : IR^ > J(T) and 


] seo = H drz eS) 
B TT) 


H d"z det (22) e- S62). 
n Ox) 


—One can do usual Monte Carlo + reweighting by regarding 


Sus r(x) := S(z(T;x)) — In L (=>) 


Ox) 


as the effective classical action. 
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Review 


Real-time dynamics 
This method is applied to Schwinger-Keldysh path integral for 


03? Ess Lä = Re[Cr LI 
T. EN — Im[Gr,odt.t)] y? 


0.0 0.5 1:0 1.5 2.0 


Feynman propagators at 8 = 0.8. Taow = 0.2. (Alexandru, et. al. (1605.08040)) 
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Summary & Conclusion 


Summary and conclusion 
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Summary and Conclusion 


e Lefschetz-thimble method gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


o One-site Hubbard model is a nice toy model to play with the 
sign problem. 


e Interference of complex phases among Lefschetz thimbles plays a 
pivotal role for the (baryon) Silver Blaze problem. 


e Recent developments may enable us to study nonperturbative 
field theories with the sign problem. 
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Backups 
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Complex Langevin and Lefschetz thimbles 


For any solutions of the DS eq, 


(0) = 2 374. | Mo, 


c o 


ind, € C. To reproduce physics, d, = (Ko, R) € Z. 


So far, we ONLY assume the convergence of the complex Langevin 
method. 
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Semiclassical limit 
Let us take h «& 1 for computing 


dE = ds [ dzeS/*O(z). 


| have NO idea how to compute the LHS. However, positivity of the 
probability density and its localization around zo's imply that 


Ieo 20 3 ODE y O. 


RHS is 


[ die SOM"O(z) > S Ole). 
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Semiclassical inconsistency 
In the semiclassical analysis, one obtains (for dominant saddle points) 


m NNNM Sieh 
Z \ Be.) l 


Co 


Co > 0. And, d, = (Ks, R) to get physics. = Inconsistent! 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 
We show that the complex Langevin is wrong if 
e There exist several dominantly contributing saddle points. 
e Those saddle points have different complex phases. 


Open question: Recently, our knowledge on the CL method is 
developing (Aarts, Seiler, and Stamatesc, 2009; Nishimura, Shimasaki, 2015, etc.). 
Can one understand our result from this viewpoint? 
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Proposal for modification 
Assume that 


_ Ko, R) 2nh Sisi) | 


S" (x. 


Because of the localization of probability distribution P, it would be 
given as 


P= GPs supp(P,)Nsupp(P;) = 0. 


Assumption means "CL — phase quenched multi-thimble approx." : 


ts = ee 


o 


2rh 
T äich 


Ee De): 
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Proposal for modification (conti.) 


If so, defining the phase function 


Lëns, H —iIm Sta lb = 
S" (Ze) 8 id Xsupp(Ps) (2,2), 


we can compute 


(Pq Sol (29). 


Ota" = 


(®(2n, 24))q 


This new one is now consistent within the semiclassical analysis. 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 

Caution: Our proposal evades inconsistency, but is not necessarily 
correct. Can we improve the proposal? 
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